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We report measurements of charmed-hadron (D", D") production cross sections at mid-rapidity 
in p + p collisions at a center-of-mass energy of 200 GeV by the STAR experiment. Charmed 
hadrons were reconstructed via the hadronic decays D° K-n+, D*+ D°tt+ R- TT^TT^ and 
their charge conjugates, covering the pr range of 0.6—2.0 GeV/c and 2.0—6.0 GeV/c for and 
D*^ , respectively. From this analysis, the charm-pair production cross section at mid-rapidity is 
da/dy\y'LQ= 170 ± 45 (stat.) Igg (sys.) fih. The extracted charm-pair cross section is compared to 
perturbative QCD calculations. The transverse momentum dierential cross section is found to be 
consistent with the upper bound of a Fixed-Order Next-to-Leading Logarithm calculation. 

PACS numbers; 25.75.-q 



I. INTRODUCTION 

The primary goal of ultra-relativistic heavy-ion 
experiments at the Relativistic Heavy Ion Collider 
(RHIC) is to search for and characterize the new 
state of matter vifith partonic degrees of freedom, 
namely the Quark Gluon Plasma (QGP), predicted 
by Quantum ChromoDynamics (QCD) In high- 
energy collisions at RHIC, heavy quarks (c, b) arc ex- 
pected to be created from initial hard scatterings Q 
and the relative changes in their masses are small by 
the strong interactions with the QCD medium Q. 
Thus they carry clean information from the sys- 
tem at the early stage. The interaction between 
heavy quarks and the medium is sensitive to the 
medium dynamics, therefore heavy quarks are sug- 
gested as an "ideal" probe to quantify theproperties 
of the strongly interacting QCD matter j^-Q. Con- 
sequently, measurements of heavy-quark production 
over a wide transverse momentum (pr) region in 
proton-proton {p + p) collisions are critical to pro- 
vide a baseline for understanding the results from 
heavy-ion collisions. In particular, precise knowl- 
edge of the total charm production cross sections 
from p + p to central heavy-ion collisions is criti- 
cal to understand both open charm and charmonium 
production mechanisms in the QGP medium formed 
in central heavy- ion collisions at RHIC 0, 01 • 

In elementary particle collisions, processes involv- 
ing heavy quarks with masses much larger than the 
QCD scale (Aqcd) are, in principle, amenable to 
perturbative QCD (pQCD) calculations. For heavy- 
quark production cross sections at large momen- 
tum transfer , Fixed-Order Next-to-Leading Log- 
arithm (FONLL) pQCD calculations, where pr ^ 
TTic, are expected to work reasonably well How- 
ever, calculations of the charm cross section at low 
Pt become complicated because charm quarks can- 
not be treated massless flavor. Furthermore, 
in the low momentum transfer region there is a 
large uncertainty in the gluon density function, and 



deceased 



the strong coupling constant increases dramatically. 
Thus, perturbative QCD calculations have little pre- 
dictive power for the total charm cross section in 
high-energy hadron-hadron collisions [loj . In view of 
these theoretical issues, experimental measurements 
become necessary and in turn provide constraints 
that improve theoretical calculations. 

Measurements of inclusive charm production have 
been carried out through two main approaches: i) 
single leptons from heavy flavor semi-leptonic de- 
cays, and ii) charmed hadrons from hadronic de- 
cays. The advantages of the first method include 
an experimentally triggerable observable and rela- 
tively large decay branching ratios, thus resulting in 
relatively large statistics. However, interpretations 
of the experimental results contain ambiguities be- 
cause a) leptons are produced by various charmed 
and bottomed hadron decays, and h) heavy-flavor 
hadrons contributing to leptons at a certain pT can 
come from a wide kinematic region due to the de- 
cay smearing. The second method suffers from a 
large combinatorial background when all particles 
from the collision vertex are included, without any 
reconstruction of the secondary weak-decay vertices. 
This background is particularly large (S/B is in the 
order of 1 : 10'^) in heavy-ion collisions. 

There are many measurements of the charm pro- 
duction cross section in low energy p + p or p + 
A collisions via both semi-leptonic and hadronic de- 
cays at CERN and Fcrmilab ^ [H]. Results for 
the total charm cross sections (from measurements 
with reasonable extrapolations) are consistent with 
Next-to-Leading-Order (NLO) pQCD calculations. 
At high energies, the Collider Detector at Fermilab 
(CDF) collaboration at the Tevatron measured the 
charmed-hadron cross sections at > 5 GeV/c in p 
+ p collisions at ^/s = 1.96 TeV, and results for 
, D'^ and D*"*" mesons are consistent with the up- 
per bounds of FONLL pQCD calculations At 
RHIC energies, charm production has been studied 
mainly via semi-leptonic decay electrons from p + 
p to Au + Au collisions |14| - [l8| . The result from 
p + p collisions is also consistent with the upper 
bound of FONLL pQCD calculations at prie) > 2 
GeV/c. Measurements of the Z?" cross section by 
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the reconstruction of hadronic decays were carried 
out in d + Au collisions IJ] , but no measurement of 
the charmed-hadron production cross section in p + 
p collisions has been made at RHIC until now. 

In this paper, we report measurements from the 
STAR experiment of the charmed-hadron (D*^, D*) 
production cross section at mid-rapidity inp + p col- 
lisions at -^/i = 200 GeV. Charmed hadrons, 1?° and 
D*, were reconstructed via hadronic decays in the 
transverse momentum ranges of 0.6—2.0 GeV/c and 
2—6 GeV/c, respectively. The pt differential pro- 
duction cross sections are compared to pQCD theo- 
retical calculations, and a total charm cross section 
is extracted. 

The paper is organized as follows: Section II de- 
scribes the experimental setup, the data set, and 
the particle identification method used in this analy- 
sis. Section III explains the hadronic reconstruction 
for and D* mesons in detail. Section IV dis- 
cusses the reconstruction efficiency, acceptance, and 
trigger/vertex corrections. Details of the systematic 
uncertainties are discussed in Section V. The trans- 
verse momentum differential production cross sec- 
tion is presented in Section VI and it is compared 
with pQCD FONLL and PYTHIA ^ calculations. 
The results are summarized in Section VII. 



II. EXPERIMENTAL SETUP 
A. Detector Apparatus 

The data used in this analysis were recorded by 
the Solenoidal Tracker At RHIC (STAR) detec- 
tor [13] ■ The STAR detector is a multi-purpose 
spectrometer with large rapidity coverage. The ma- 
jor subsystems at mid-rapidity sit inside a solenoidal 
magnet which provides a uniform magnetic field of 
0.5 Tesla along the beam axis. Subsystems used 
in this analysis are the Time Projection Chamber 
(TPC) the Time-Of-Fhght (TOF) detector [H, 
the barrel and endcap Electromagnetic Calorimeters 
(EMC) [H,!!!] , and two trigger detector subsystems: 
the Vertex Position Detector (VPD) ^ and the 
Beam Beam Counters (BBC) |26l |. 

The TPC is the main tracking detector, covering 
the full azimuthal angle at pseudo-rapidity jryj < 1 
for tracks crossing all 45 padrows [2l|. It measures 
the chargcd-particle momenta and provides particle- 
identification (PID) capability via the ionization en- 
ergy loss (dE/dx) in the TPC gas, allowing a clean 
separation between charged kaons and pions up to 
momentum p ^0.6 GeV/c. The barrel TOF detector 
is a newly installed subsystem, utilizing the multi- 
gap resistive plate chamber technology [l^. The 
full system consists of 120 trays covering the full az- 



imuth at |?7| < 0.9 surrounding the TPC cylinder. 
In the year 2009 run, 84 trays out of 120 for the 
full barrel were installed and used for this analysis. 
The TOF detector uses the timing recorded in the 
forward VPD detector as the start time to calculate 
the particle time of flight, which is combined with 
the momentum from the TPC to identify particles. 
The timing resolution of the TOF system, includ- 
ing the start timing resolution in = 200 GeV p 
+ p collisions, is about 110 ps, allowing separation 
of K and tt up to p ~1.5 GeV/c. The barrel and 
endcap EMCs are designed to identify electrons and 
photons, covering the full azimuthal angle at \r]\ < 
1 and 1 < 77 < 2, respectively (23l . [23 |. They are 
fast-response detectors (< 100 ns), and were used to 
suppress the TPC pilcup-track contribution in the 
event- vertex finder by matching with charged tracks 
from the TPC. 

In addition to providing the start time for the bar- 
rel TOF detector, the VPD detector is also one of 
the trigger detectors in STAR. It has two parts sur- 
rounding the beam pipe, located on the east and 
west sides, 5.7 m away from the center of the STAR 
detector and covering 4.24 < jryj < 5.1 [25|. The 
minimum-bias (MB) trigger was defined as a coinci- 
dence signal in the east and west VPD detectors and 
a selection was made on the vertex position along the 
beam axis (Vz) to be within 40 cm of the center of 
the STAR detector. The BBC ^ consists of two 
identical counters located on each side of the TPC 
covering full azimuth and 2.1 < \ri\ < 5.0 in pseudo- 
rapidity. Each part consists of a set of hexagonal 
scintillator tiles grouped into a ring and mounted 
around the beam pipe at a distance of 3.7 m from 
the center of STAR. The BBC detector had been 
used to define the main minimum-bias trigger in p 
+ p collisions before the minimum-bias trigger was 
used in 2009. A small sample of BBC minimum-bias- 
triggered events were collected in 2009 to check for 
a trigger bias. Details of the minimum-bias trigger 
bias and correction will be discussed in Sect. IV. 



B. Data Sets and Event Selection 

The data sample used in this analysis consisted of 
minimum-bias-triggered p + p collisions at y/s ~ 200 
GeV, recorded in 2009 by the STAR experiment at 
RHIC. 

The intrinsic drift time for electrons from the 
center to one end of the TPC is on the order of 
40 fis. Thus, in high-luminosity p + p collisions, 
one TPC event usually contains tracks from col- 
lisions originating from non-triggered bunch cross- 
ings. These "pileup events" will lead to additional 
tracks recorded in the TPC, in addition to those 
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from the triggered event. This effect was not sig- 
nificant in previous RHIC runs, but the increase 
in the colhsion rate during 2009 to several hundred 
kHz made this a significant effect. The T4 position 
from offline VPD data has a resolution of 2.5 cm 
for minimum-bias events, which can provide a use- 
ful constraint to select the real event that fired the 
trigger. Figure [U upper panel, shows the correla- 
tion between the Vz positions from the TPC and 
the VPD. Events with TPC vertices along the di- 
agonal correlated band are real ones that fired the 
VPD minimum bias trigger. In Fig.[TJ bottom panel, 
the solid black histogram shows the 1-D differ- 
ence between the first TPC-determined vertex posi- 
tion and VPD-determined vertex position. By ap- 
plying a Vz difference cut | AT^zj < 6 cm, most of the 
TPC pileup events can be removed. There still re- 
main random associated correlations that enter into 
this cut window (~7% level, calculated using a two- 
Gaussian fit). To further suppress this contamina- 
tion, we required the TPC event vertices to have at 
least two tracks that match with hits in the barrel 
and endcap EMCs (this vertex is treated as a "good" 
vertex). The red dashed histogram in Fig.[TJ bottom 
panel, shows the AVz distribution after this selec- 
tion. The random associated pileup events in the Vz 
difference cut window arc now suppressed to ^^2% of 
the total, while the corresponding loss of real events 
is ^--^15%. In total, 105 million minimum-bias events 
were used in the charmed-hadron analysis. 



C. Track Reconstruction and Particle 
Identification 




-40 -20 20 



A V, = vr'=-Vr° (cm) 

FIG. 1: (color online) Upper Panel: Correlation of V^^^ 
versus VY^^- Bottom Panel: AVz distributions. A 
"good" vertex requirement rejects most of the pileup 
events. Blue vertical lines indicate the cuts for the 14 
selection. 



Charged particle tracks are required to point 
within \ri\ < 1 in order minimize TPC acceptance 
effects during reconstruction. Tracks must have 15 
out of a maximum of 45 points used in track fitting 
(nFitPts), and at least 52% of the total possible fit 
points in order to avoid double-counting split tracks. 
Tracks are required to have a distance-of-closest- 
approach (DCA) to the collision vertex of less than 2 
cm to suppress background tracks produced by sec- 
ondary scattering in the detector and also long-lived 
particle decays. The STAR track pointing resolu- 
tion with the TPC alone does not have the precision 
to separate charm secondary decay vertices from the 
collision vertices. 

Particle identification for final-state charged 
hadrons was carried out with a combination of 
dE/dx in the TPC and the particle velocity {(3) mea- 
surement from the barrel TOF detector. Thus the 
normalized dE/dx [na'^^'^^) and 1//3 (ncr™^) dis- 
tributions were used to select daughter particle can- 
didates. They are defined as follows: 



dE/dx 



. UE/dxTZ 

dE/dx^^ 

RdE/dx 
1 ]_ 



//3 



(1) 

(2) 



where the superscripts "mea" and "th" are mea- 
sured and theoretical values, respectively. The X 
denotes expected values which are calculated with 
respect to one kind of particle species (tt or K). 
RdE/dx and Ri/[j are the experimental dE/dx and 
1//3 resolutions, respectively. With the above defini- 
tions, the two resulting distributions can be approx- 
imated by Gaussian distributions with mean^O and 
(7 ~1). Figure [2] shows the na 
na 



dE/dx 

™^ distributions versus particle momentum. 



dE/dx 
K 



and 



Daughter kaon (pion) candidates are selected by 
requiring |n(T^^'^''^| < 2 (\na'i^^'^^\ < 2). In addi- 
tion, to improve the significance of the reconstructed 
D'^ signal, the kaon daughter tracks were required to 
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1 2 3 4 5 

Momentum (GeV/c) 



dE/dx 



FIG. 2: (color online) Distributions of not^^'^^ , ^'^'k 
and najf^ versus momentum are shown in panels (a), 
(b), and (c), respectively. The latter is shown after 
dE/dx cuts were applied. 



have a valid hit in the TOF detector and then se- 
lected with a TOF PID cut, which is denoted as the 
red dashed lines in panel (c) of Fig. [21 In order to 
have good efficiency and considering pion identifica- 
tion is good enough with dE/dx only, we did not 
require pion to match with TOF. 



III. CHARMED-HADRON 
RECONSTRUCTION AND RAW YIELD 
EXTRACTION 

A. D'^ Reconstruction 

_D" and mesons were reconstructed via the 
hadronic decay D^{D'^) — > K^n^ with a branching 
ratio of 3.89%. The analysis technique is the same as 
that used for a analysis in -I- Au collisions [3] . 
Inp + p collisions, the mixed-events technique is not 
suitable for describing the background due to large 
contribution of correlated jets. Therefore, two dif- 
ferent techniques were used to reproduce the back- 
ground: the like-sign and track-rotation methods. 
Since the tt~ and 7r+ production is symmetric in 
the STAR uniform acceptance and their yield ra- 



tio is measured to be 0.988 ± 0.043 l27|, the like- 
sign (LS) method is used and a pair combination 
with the same charged sign is expected to repro- 
duce the background without the signal correlation. 
The opposite-sign backgrounds, which go into the 
residual background, arc only several percent of total 
background and will be discussed later. The track- 
rotation (Rot) technique has been used in many 
measurements [2^ . This method is based on the as- 
sumption that by rotating the daughter kaon track 
by 180 degrees in azimuth, the decay kinematics are 
destroyed. Thus the invariant mass distribution af- 
ter rotation is able to reproduce the random com- 
binatorial background. Figure |3] shows the invari- 
ant mass distributions of Kit candidates. Panel (a) 
shows the invariant mass distributions for Kt: pairs 
(0.6 < priKn) < 2.0 GeV/c) with unlike-sign (US) 
before background subtraction, with like-sign, and 
with rotated kaon momentum. The distributions 
from the like-sign and track-rotation techniques de- 
scribe the background well. Panel (b) is the unlike- 
sign Kit invariant mass distribution after combinato- 
rial background subtraction. A significant i^*(892) 
peak is observed. The secondary small peak at about 
1.4 GeV/c^ is the i4r2(1430). A direct zoom-in view 
of the vicinity around the mass region is shown 
in Fig. U (panel (a) for subtraction of like-sign back- 
ground, and panel (b) for the rotational case). Solid 
symbols depict the same distributions as shown in 
Fig. [3] and Fig. [S] in two different Z?" pr bins. One 
can see there is still some "residual" background af- 
ter like-sign or rotational background subtraction. 
The possible sources to the residual background have 
been investigated using PYTHIA simulations. We 
performed the same reconstruction as we did on 
the data, for the foreground and background distri- 
butions. From these simulations, we have learned 
that the possible sources that can contribute to 
this residual correlated background include: corre- 
lated hadron pairs from decays (mostly resonances) 
where the real daughters were mis-identified as Kn 
pairs; Ktt pair from other decay channels of D'^ (e.g. 
_ft'~7r+7r°) where the other daughters are missed in 
the reconstruction; same-charge K~t:~ pairs from 
multi-body decays of D° — >■ K^'!t^tt^tt~; Ktt pairs 
from jet fragmentations; etc. The different shape 
of the residual background from LS and Rot back- 
ground subtraction in the data can be qualitatively 
reproduced by PYTHIA simulation. The magni- 
tude of the residual background depends on how to 
choose the normalization for the like-sign or rota- 
tional background, as qualitatively understood from 
the PYTHIA simulations. However, the change of 
the residual background magnitude due to different 
normalizations has a very small impact on the final 
extracted signal counts, and it has been included in 
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the systematic uncertainties. We used an empirical 
polynomial function to describe it and the choice of 
this empirical function was also included as one of 
the systematic source to the raw yields. A Gaussian 
function is used to fit the signal. The raw yield of the 
is obtained by fitting the data (blue solid circles) 
with a fit function representing the sum of signal and 
background (red dashed curve) in the mass region of 
1.72 < Mk^ < 2.05 GeV/c^. The signal after the 
residual background subtraction is shown as the red 
open circles. The Gaussian function used to describe 
the signal is shown as the blue dashed curve. The 
total signal consists of 4085 ± 938 counts. 

The signals after background subtraction for two 
Pt bins are shown in Fig. [5] Panels (a), (c) and (b), 
(d) show the signals from LS and Rot background 
subtraction, respectively. The raw yields and 
statistical errors extracted from the two background 
methods arc listed in Table ID The average values 
of the counts from the LS and Rot background 
methods are used to calculate the final raw yield 
in each pT bin. The mean and width from the Gaus- 
sian fits are compared with MC simulation in Fig. [S] 
(left panels). The single D'^ and D* are embedded 
into the real data and simulated in the full STAR 
GEANT reconstruction chain, taking into account 
detector response and material effect. The sig- 
nal mean value from an open-parameter fit shifts to 
lower mass due to kaon energy loss at low p^, which 
is not fully accounted in the simulation due to possi- 
bly missing material budget. The systematic uncer- 
tainty in determining the raw yields as well as 
the potential double-counting issue due to particle 
misidentification will be discussed in Sect. V A. 



TABLE 1: D° raw yields. 



Pt range (GeV/c) 


0.6-1.2 


1.2-2 


PT (GeV/c) 


0.908 


1.57 


raw yields x 10^ (Rot) 


2.45 ± 0.66 


1.65 ± 0.63 


raw yields x 10^ (LS) 


1.67 ± 0.74 


2.40 ± 0.64 



B. D* Reconstruction 

D*^ mesons were reconstructed via the decay se- 
quence D*+ D°7r+ {BR = 67.7%), D° ^ K-Tr+ 
and its charge conjugate. We followed the same 
analysis technique as described in Ref. [2§|. The 
daughter particles were still identified by dE/dx in 
the TPC because a) most of the D* decay daugh- 
ter particles that fall inside the STAR acceptance 
with higher momenta are located in the region where 
the TOF PID improvement is very limited; and b) 
the signal suffers significant losses due to incomplete 




M(Kjt) (GeV/c' 



M(Kjt) (GeV/cT 



FIG. 3; (color online) (a) Invariant mass distributions of 
raw Kn combinations for unlike-sign pairs (circles) , like- 
sign pairs (triangles), and kaon momentum rotated pairs 
(line), (b) Residual distributions after subtracting the 
like-sign distribution (triangles) and rotation pair distri- 
bution (dots) from the unlike-sign distribution. 
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FIG. 4: (color online) Invariant Kir mass distributions in 
the _D° mass region after like-sign (a) and track-rotation 
(b) background subtraction. Solid circles show the sig- 
nal and a residual background. A Gaussian function and 
a 2nd order polynomial function were used to describe 
the signal and residual background, respectively. Open 
circles show the signal after residual background subtrac- 
tion. 



TOF acceptance in 2009. Compared to the cuts used 
in Ref. [1^ , the pr threshold cut for the 7r+ (from D* 
decays), denoted as 7r+, was lowered to 0.15 GeV/c. 
The ratio, r, of transverse momenta from the D° 
and 7r+ was required to be 7 < r < 20. These two 
changes were implemented to improve the statistics 
near the lower bound in pT- The remainder of the 
analysis cuts were the same as those used in Ref. [2§| . 

The invariant mass difference AM = M{KTnr) — 
M{Ktt) was calculated in reconstructing the D* sig- 
nal to take advantage of the partial cancellation in 
the detector resolution in measured mass distribu- 
tions. The AM distributions are shown in the up- 
per panel of Fig. [T] The "right-sign" combinations 
K^TT^nf were used to select the D*"^ candidates. 
Two independent methods - "wrong-sign" combina- 
tions i^'^TT^Tr^ and D° "side-band" combinations 
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FIG. 5: (color online) Raw signals in different pr 
bins after like-sign (a)(c) and track-rotation (b)(d) sub- 
traction. 
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FIG. 6: The mean and width from Gaussian fit to data 
(symbols) compared with MC simulations (bands) for 
and D* are shown in left and right panels, respectively. 



area indicates candidate mass selection in the 
Kmr right-sign and wrong-sign combination recon- 
struction. The line-hatched area indicates the 
side-band region (1.72 < M{KTT)/{GeY/c^) < 1.80 
or 1.92 < M{KTr)/{GcV/c'^) < 2.00) used in side- 
band combinatorial background reconstruction for 
D* . The side-band combinatorial background was 
used to obtain the raw D* yields for better statis- 
tics and also because side-band distributions do not 
suffer from the double-counting issue due to par- 
ticle misidentification. The difference between the 
yields obtained from the side-band method and the 
"wrong-sign" method was included in the system- 
atic uncertainties. Details in determining the uncer- 
tainties on the raw D* yields including the double- 
counting effect will be discussed in Sect. V A. The 
D* raw yields are summarized in Table HIl 



TABLE II: D* raw yields. 



PT range (GeV/c) 


2-3 


3-4 


4-5 


5-6 


PT (GeV/c) 


2.45 


3.44 


4.45 


5.45 


raw yields 


209 ± 58 


98 ± 35 


27 ± 11 


12.3 ± 4.1 



To obtain the cross section, the event-selection 
criteria described in the previous section were ap- 
plied. The raw distributions were further divided 
into pt slices to obtain the raw D* yields in each 
bin. Figure [5] shows the D* candidates and back- 
ground distributions in different pT bins. The bot- 
tom panel on each plot was generated by subtracting 
the "side-band" background from the "right-sign" 
candidates. The mean and width from Gaussian fits 
arc compared with MC simulation in the right panel 
of Fig. m and it shows the obtained D* peak posi- 
tions and widths agree with the MC simulation well. 
From this analysis, the total signal consisted of 364 
± 68 counts, and the raw yield ratio of D*^ / D*^ is 
0.93 ± 0.37. 



IV. EFFICIENCY AND TRIGGER/VERTEX 
BIAS CORRECTION 

The final charmed-hadron cross section in p + p 
collisions is calculated as follows: 



- were used for combinatorial background recon- 
struction. The plot illustrates that both methods 
reproduce the combinatorial background very well. 
The events displayed in this figure are all minimum- 
bias events without event vertex selections, which 
demonstrates the significance of D* signal. The 
lower panel in Fig. [7] shows the Kn invariant mass 
distribution after requiring the D* candidate cut 
(0.144 < AM < 0.147 GeV/c^). The cross-hatched 



E— = — — — gNSD r 

dp3 2tt ■ e^ec ' BR ■ PT^prAy ' Nmb ' 

(3) 

where unsd is the total Non-Singly Diffractive 
(NSD) cross section, which is measured at STAR to 
be 30.0 ± 2.4 mb [s^- A^mb is the total number of 
minimum-bias events used for the analysis. ANjj is 
the raw charmed-hadron signal in each pT bin within 
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FIG. 7: (color online) Upper: Raw D* candidate signal 
from the right-sign combinations in all p + p minimum- 
bias events. Histograms are combinatorial background 
distributions from "wrong-sign" and "side-band" meth- 
ods. Lower: Raw D" candidates after requiring the D* 
candidate cut (0.144 < AM < 0.147 GeV/c^). 



FIG. 8: (color online) Raw D* signals in different pr 
bins. In each plot, the bottom panel distribution is gen- 
erated by subtracting the side-band background from the 
right-sign distribution. Variable binning is used in the 
bottom panel for better illustration. 



a rapidity window Ay. BR is the hadronic decay 
branching ratio for the channel of interest. There 
are two correction factors: Crec, which is the recon- 
struction efficiency including geometric acceptance, 
track selection efficiency, PID efficiency, and analy- 
sis cut efficiency; and fug.vtxipr), which is the cor- 
rection factor to account for the bias between the 
minimum-bias sample used in this analysis and the 
total NSD sample. This bias is mainly caused by the 
VPD trigger and event vertex reconstruction, and it 
may have a dependence on the charmed-hadron pr- 
In the following sections of the paper, the condi- 
tion that requires the event to fire the VPD trigger 
and to have a good vertex will be referred to as the 
"Analysis Condition" . 



A. Reconstruction Efficiency 

The reconstruction efficiency for charmed hadrons 
was obtained by embedding Monte Carlo (MC) 



simulated charmed-hadron tracks into the real 
minimum-bias events. The MC charmed-hadron 
tracks were processed through a full GEANT detec- 
tor simulation [3l| with a representation of the 2009 
STAR geometry. The raw detector-response signals 
were mixed together with those from the real data 
and processed through the full STAR offline recon- 
struction chain to obtain the detector response ef- 
ficiency in a realistic environment. The input MC 
track multiplicity was constrained to have negligible 
effect on the final tracking efficiency due to increased 
occupancy in the TPC. 

Figures 1^ and [TUl show the D'^ and D* reconstruc- 
tion efficiency versus within |y| <1. In Fig.|9l the 
solid squares denote the reconstruction efficiency for 
both daughters selected and identified by the TPC, 
while the solid circles denote the reconstruction effi- 
ciency with additional PID selection from the TOF 
detector for the kaon daughter. The combined TOF 
efficiency, including the acceptance, matching be- 
tween TPC tracks and TOF hits, and PID selection 
efficiency, is around 45% studied from the data in 
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FIG. 10: Total D* reconstruction efficiency versus D* 
Pt- 



VPD trigger study, we first compared the VPD trig- 
ger efficiencies (from tlie BBC triggered minimum- 
bias sample) from MC simulation and real data. The 
BBC trigger has been well studied and was used to 
calculate the p + p NSD cross section [l6|. Figure [TT] 
shows the comparison of the VPD trigger efficiency, 
with the requirement that there is a BBC trigger 
and a good vertex. The efficiency is studied as a 
function of the charged hadron px- The real data 
used are BBC triggered minimum-bias events taken 
in 2009 during a very low luminosity run, which min- 
imizes TPC pilcup tracks. Figure [TT] shows that 
the efficiency goes down with increasing px of mid- 
rapidity particles indicating an anti-correlation be- 
tween mid-rapidity particle production and forward 
VPD triggering. Most importantly, within the mo- 
mentum range under study, the PYTHIA MC sim- 
ulation agrees well with the data. This agreement 
provides confidence in using PYTHIA simulations to 
evaluate this correction. 
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FIG. 11: VPD trigger efficiency comparison between 
Data and Monte Carlo versus charged particle pr in BBC 
MB conditions. 



B. Trigger and Vertex Bias Corrections 

The trigger and vertex bias corrections were stud- 
ied by simulating PYTHIA events [l^ processed 
through the full GEANT detector response and of- 
fline reconstruction. The PYTHIA generator ver- 
sions 6.205 and 6.416 were both used in this study. 
We chose the PYTHIA version 6.205 with minimum- 
bias processes selected and with the CDF TuneA set- 
tings [32[ to give the centroid value of the correction 
factor because it gives better description for the par- 
ticle production in the forward rapidities than the 
6.416 version [s^. The differences between the two 
versions as well as different parameter settings have 
been included to estimate the systematic uncertainty 
of the trigger and vertex bias correction factor. 

To validate the PYTHIA generator in simulat- 
ing particle production in the forward region for the 



The correction factor /trg,vtx can be related to the 
ratio (iVu/iVinb) for the pure minimum-bias condi- 
tion and the "Analysis Condition", i.e. 



/trg,vtx(PT) 



Nd{pt)/N, 



mb 



,rtrg,vtx 



tig,vtx ■ 



(4) 



Two simulation samples were generated to ob- 
tain the correction factor. One sample consisted of 
PYTHIA-simulated p + p events and was used to 
obtain the fraction of minimum-bias events that sat- 
isfy the "Analysis Condition" : N^^^'^^"" /N,^h- This 
fraction was found to be 12.7% from this PYTHIA 
simulation. The other simulation sample was gen- 
erated using the same PYTHIA settings, but only 
events with at least one charmed hadron were saved 
to enhance the statistics. This sample was used to 
obtain the fraction of charmed-hadron signals that 
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satisfy the "Analysis Condition" - TV^'s.vtx^^^^ 
also studied this fraction as a function of charnied- 
hadron pT- Figure [T^] shows the calculated effi- 
ciencies for D* from different event-selection crite- 
ria. The BBC coincidence study provides a baseline 
for this simulation, which demonstrates consistency 
with previous STAR results [lOl- As expected, the 
vertex finding efficiency increases with increasing ■ 
The VPD trigger efficiency shows an anti-correlation 
with increasing D* pT, similar to that observed with 
increasing charged-hadron px- The final efficiency 
(with requirements for both vertexing and VPD trig- 
gering) is almost fiat versus leveling off at ~19%. 
The simulation for hadrons shows very similar re- 
sults. Figure [T51 shows the correction factor, /trg,vtx, 
for cross section calculations for and D* . 
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FIG. 12: D* efficiency versus D* pr with different event 
selection criteria. 
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hadron pr for cross section calculations for and D* 



V. SYSTEMATIC UNCERTAINTIES 

Sources that contribute to the systematic uncer- 
tainties in the final D-meson cross sections include: 
a) uncertainty in determining the raw ZJ-meson 
yields; b) uncertainty in determining the reconstruc- 
tion efficiency; c) uncertainty of the total NSD cross 
section and d) uncertainty in determining the trig- 
ger/vertex correction factor. Uncertainties due to 
particle identifications will enter in both a) and b) 
which will be discussed in the following subsections. 
We consider a) as point-by-point uncorrelated sys- 
tematic uncertainties. Although b) is correlated in 
Pt, it is not simply a normalization uncertainty, and 
the exact correlation in px is not known. Therefore 
we include b) in the point-by-point uncorrelated sys- 
tematic uncertainties. Finally, c) and d) are overall 
normalization uncertainties. 



A. Uncertainty in Raw Yields 

Different choices on background reconstruction 
methods, function fits and mass binning were used 
to evaluate the systematic uncertainty in the raw 
ZJ-meson yields. In the D° analysis, the difference 
between the yields extracted from Rot and LS meth- 
ods is 15.6-18.9%. Fitting the D° peak with fixed 
parameters from simulation estimates lower yields 
of 28.2% and 6.1% for the two D° pr bins. The 
systematic uncertainties from different mass bin- 
ning and different fit regions are estimated to be 
'^5—7%. The systematic uncertainties in determin- 
ing the raw D* yields include contributions from the 
difference obtained between the "side-band" and the 
"wrong-sign" methods, and the difference between 
bin counting and Gaussian fitting methods, varying 
^6—11% in the pr range 2—6 GeV/c. The choice 
of mass binning and fitting range had a negligible 
effect on the extracted yields. 

In D'^ meson reconstruction, if the kaon (pion) 
daughter is misidentified as a pion (kaon) then two 
daughters from a real decay will show up as ad- 
ditional combinations with a wider mass distri- 
bution due to wrong mass assignments. Thus one 
signal will be counted twice; once as a and 
again as a Z?". A Monte Carlo simulation was used 
to evaluate the fraction of such double counting oc- 
currences in the reconstruction. Based on re- 
alistic dE/dx and TOF FID resolutions extracted 
from real data, the probability that kaons (pious) 
can be misidentified as pions (kaons) at a given pt-, 
using these FID selections, was obtained. Assuming 
a D*^ candidate, this procedure provides an estimate 
of the probability that both daughters are misiden- 
tified and then reconstructed as a In Fig. [Ml 
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the open and closed circles show the double-counting 
fraction, relative to the total real signal, for two dif- 
ferent PID selections: a) both daughters are iden- 
tified by TPC dE/dx; b) the kaon daughters arc 
identified by the TOF, while pious arc identified by 
the TPC. The sharp increase at very low pT (iden- 
tifying both daughters using dE/dx) is due to the 
case where a decays almost at rest [px ^0), and 
the two daughters are produced in the momentum 
region where the kaon and pion dE/dx bands cross, 
therefore maximizing the misidentification probabil- 
ity. The plot shows that when the kaon daughter is 
identified by the TOF, the double-counting fraction 
is negligible in our pr coverage region (0.6—2.0 
GeV/c). 
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FIG. 14: D double-counting fraction due to particle 
misidentification in two PID selections and D* wrong- 
sign over counting fraction versus D-meson pr- 

Double counting the D'^ may also impact recon- 
struction of D*. However, the impact is different 
because of a charge sign requirement on the soft pi- 
ous. If both daughters from a D° are misidentified 
(D° is reconstructed as D^), then the combination 
from the same signal will become K~^Tr~T:'^. It will 
not contribute to the right-sign distributions, but in- 
stead, will enter into the "wrong-sign" (background) 
distributions if the mass also falls into the (^°) 
mass selection window. Thus the double counting 
in "wrong-sign" background will contribute to an 
undercounting in the total signal if the wrong-sign 
background is subtracted from the right-sign distri- 
bution. Since the right-sign combination was also re- 
quired, the misidentification does not affect the side- 
band background distributions. In the real analysis, 
the side-band background subtraction was used to 
extract the raw signal, but also the difference be- 
tween side-band and wrong-sign methods was used 
for systematic uncertainty estimation. Since the 
wrong-sign distribution can be overestimated due to 
particle misidentification, the systematic error from 



the difference between the two methods would be 
overestimated. This was avoided with better under- 
standing of the wrong-sign overcounting. The red 
triangles in Fig. [U denote the over counting frac- 
tion in the D* wrong-sign background to real sig- 
nals. It is very close to the D° double-counting frac- 
tion, since they are from the same source. The slight 
difference comes from the additional candidate 
selection cuts used in the D* reconstruction. This 
fraction was used to compensate for the difference 
between the two background methods and as a way 
to improve the assessment of the systematic uncer- 
tainties in the extraction of the raw D* yields. 



B. Uncertainty in Reconstruction Efficiency 

The systematic uncertainties of the reconstruction 
efficiencies were obtained following similar methods 
used in other particle cross section measurements by 
changing the daughter track selection criteria and 
comparing the difference between the data and the 
MC. In this analysis, it was studied by changing the 
minimum number of fit points (nFitPts) in the TPC 
from 15 to 25 and the DCA to the collision vertex 
from 2 cm to 1 cm. The uncertainty was then quan- 
tified by the difference in the remaining fractions 
after cut changes between the data and the MC. For 
each cut change, the uncertainties were calculated 
for each decay daughter and added together linearly 
to obtain the total for and D* . The systematic 
uncertainties on the PID cut efficiencies (from both 
dE/dx and TOF) were estimated to be <1% and 
neglected in the total uncertainty. Then the uncer- 
tainties from the cut changes on uFitPts and DCA 
were added in quadrature to obtain the total sys- 
tematic uncertainty on the reconstruction efficiency. 

The point-by-point systematic errors including 
uncertainties in raw yields and reconstruction effi- 
ciency for the -D" and D* cross sections in each pr 
bin are summarized in Table IIIII 



TABLE III: D° (0.6-2 GeV/c) and D' (2-6 GeV/c) 
point-by-point systematic errors (%) 



PT (GeV/c) 


0.6-1.2 


1.2-2 


2-3 


3-4 


4-5 


5-6 


raw yields 


4-18.9 
-33.9 


+15.6 
-16.8 


9.4 


6.5 


11.0 


6.6 


nPitPts 15->25 


3.8 


3.2 


7.2 


4.7 


5.9 


4.7 


DCA 2->l (cm) 


6.6 


7.1 


13.6 


12.7 


11.6 


10.7 


quadratic sum 


+20.8 
-34.8 


+17.8 
-18.5 


18.1 


15.1 


17.1 


13.5 
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C. Overal Normalization Uncertainty 



The overall normalization uncertainty for the to- 
tal NSD cross section has been studied before and 
reported in previous STAR publication [sO]- It 
was estimated to be 8.1%, including the uncer- 
tainty from measuring the absolute BBC cross sec- 
tion and that of BBC triggering efficiency. The un- 
certainty from the trigger/ vertex bias correction fac- 
tor amounts to 5.2% by varying different PYTHIA 
versions (6.205 vs. 6.416) and different parameter 
settings in the simulation. We also considered the 
impact from pileup TPC tracks as an additional sys- 
tematic source on the correction factor , and the un- 
certainty was estimated to be 4.0% by comparing 
the result with a conservative luminosity level for 
this data set to that from pure PYTHIA simulation 
without pileup. 

These uncertainties were added in quadrature, 
which gives 10.4% overall normalization uncertainty 
for the D-meson cross sections. 



VI. RESULT AND DISCUSSION 

After the reconstruction efficiency and trig- 
ger/vertex bias correction factor were applied, the 
differential production cross sections for D° and D* 
m p + p collisions at ^/s = 200 GcV were ex- 
tracted, as shown in Fig. [151 The vertical bars 
on the data points indicate the statistical uncer- 
tainties, while the brackets indicate the bin-to-bin 
systematic uncertainties described in previous sec- 
tion. The _D" and D* cross sections were divided by 
the charm quark fragmentation ratios 0.565 ± 0.032 
(c ^ and 0.224 ± 0.028 (c ^ respec- 
tively, to convert to the cc production cross section. 
The charm quark fragmentation ratios are measured 
from CLEO and BELLE experiments near the T res- 
onance 



3J] . The uncertainties of the fragmentation 



ratios are taken into account as systematic errors 
in calculating the cc production cross section. A 
power-law fit to the data points was performed with 
the following function [l4| : 



^d?a _ da 2{n - l)(n - 2) 
^p3 j„. o^s>/„ \2 



Pt 



{pT){n-i)/2 



)-". 
(5) 

and shown as the solid red line in the figure. The 
fit quality with the power-law function, measured 
as x^/ndf, is 0.9/3 with statistical errors and 3.7/3 
with point-by-point systematic errors, respectively. 
The latter was used to extract the systematic un- 
certainty on the Pt integrated cross section from 
point-by-point systematic sources. The obtained cc 
production cross section at mid-rapidity is. 



^ If^o = 170 ± 45 (Stat.) (sys.) ^b. (6) 

The term with sys. includes the uncertainty aris- 
ing from the bin-to-bin systematic uncertainties and 
from the extrapolation to the low-pT region, which 
is not measured. The FONLL upper limit and 
PYTHIA-l-tune fits are used for the \avj-pT extrapo- 
lation, which gives -1-6.2% and -16.4% uncertainties, 
respectively. At mid-rapidity, about 67% of the D 
meson yield falls in the measured pt region. The 
mean transverse momentum of charmed mesons is 
found to be 1.06 ± 0.14 (stat.) ± 0.09 (sys.) GeV/c. 
The charm-pair cross section at mid-rapidity from 
this measurement is consistent with STAR' s pre- 



vious measurement in d-l-Au collisions |14l | at 1.7cr 
(cr is the averaged total uncertainty between two re- 
sults), providing negligible nuclear effects in d-l-Au 
collisions. 

Also shown in Fig. [15] are the upper and lower 
edges (blue dashed lines) of a FONLL pQCD cal- 
culation taken from Ref. [1]. Our results are con- 
sistent with the upper limit of the FONLL pQCD 
calculation in a wide pT region. It is observed 
that the charmed-hadron cross sections measured 
by CDF and ALICE [H at energies up to 7 
TeV are also close to the upper limits of FONLL 
pQCD calculations. This may help set constraints 
on the parameters used in the FONLL calculations, 
e.g. on the choice of renormalization or factoriza- 
tion scales, which are the main parameters varied to 
obtain the upper and lower limits on these calcula- 
tions. However one should note the valid pr region 
of FONLL calculations when applying such a analy- 
sis since FONLL calculations are supposed to work 
when Pt S> mc- 

The charm cross section at mid-rapidity was ex- 
trapolated to full phase space using the same extrap- 
olation factor, 4.7 ± 0.7, as in a previous publication 
[3], and the extracted charm total cross section at 
= 200 GeV is 

cTcc = 797 ± 210 (stat.)i208 (gys.) ph. (7) 

Shown in Fig. I16[ the data were also compared 
with PYTHIA calculations. PYTHIA version 6.416 
was used as it has been tuned to describe the mid- 
rapidity Tevatron data. We tried PYTHIA calcula- 
tions with the following sets of parameters to com- 
pare with our measurements: 

a) Default MSEL = 1. 

b) PHENIX tunc: MSEL = with MSUB(11,12, 
13, 28, 53, 68) on, PARP(91) ((fc^)) = 1.5 
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FIG. 15: (color online) cc production cross section as 
inferred from D*^ and D* production in p + p collisions at 
-^5 = 200 GeV compared with FONLL calculations. The 
D° and D* data points were divided by the charm quark 
fragmentation ratios 0.565 (c — >• D°) and 0.224 (c — >■ 
D*^) [sj, respectively, to convert to the cc production 
cross section. 



FIG. 16: (color onhne) cc production cross section as 
inferred from and D* production in p + p coUisions 
at ^/s = 200 GeV compared with PYTHIA calculations. 
Data are fitted with PYTHIA spectra with a overall scale 
parameter for the purpose of shape comparison only. 



VII. SUMMARY 



GeV/c, MSTP(32) {Q^ scale) = 4, CKIN(3) 
(min. partonpj_) = 2 GeV. 

c) This tune: MSEL = 1, PARP(91) ((fc_L)) = 
1.0 GeV/c, PARP(67) (parton shower level) = 
1.0. 

The choice of modifying the primordial {k±) (the 
Gaussian width of primordial kr in hadrons) and the 
parton shower level parameters from default values 
(2 GcV/c and 4, respectively) in this tunc was sug- 
gested by the matching of scales in heavy-flavor pro- 
duction at lower energies [36j , which has been noted 
in PYTHIA [H]. The CDF tune A parameters 
which were tuned to reproduce mid-rapidity jet and 
"underlying event" results at Tevatron energies, are 
included as defaults in PYTHIA v6.416. "PHENIX 
tune" parameters are those used in the PHENIX 
charm continuum contribution estimation from di- 
electron measurements [s^l • The default parton dis- 
tribution function (CTEQ5L) was used in all three 
cases. 

All ground-state charmed hadrons (1?°, 
and A+) were added together in the rapidity window 
\y\ < 1 to obtain charm cross sections. The data 
was then fitted with the PYTHIA calculations with 
a overall scale factor as the unique free parameter. 
The charm production pT spectrum with this tunc 
gives best x^- 1-41 (this tunc), 4.97 (dcfauh), 5.96 
(PHENIX tune). This is the first direct D- meson 
measurement that goes down to such a low px, which 
constrains the model parameters better. 



In summary, measurement on the charmed me- 
son (I?" and D*) production cross sections via their 
hadronic decays in p + p collisions at \/s = 200 
GeV has been reported. The charm pair produc- 
tion cross section at mid-rapidity extracted from 
this analysis is da/dy\f^Q^ 170 ± 45 (stat.) tfg 
(sys.) /xb. The charm total cross section at ^/s 
= 200 GeV is estimated as 797 ± 210 (stat.) tf^l 
(sys.) /ib. The reconstructed charmed mesons 
cover the pr range from 0.6-6 GeV/c. The charm- 
pair transverse momentum differential cross sections 
from this analysis are consistent with the upper 
bound of a Fixed-Order Next-to-Lcading Logarithm 
perturbativc QCD calculation. When comparing to 
PYTHIA model calculations, we found that a cal- 
culation with smaller primordial {k±) and parton 
shower level compared to CDF TuneA settings de- 
scribes the shape of the pr distribution of data. 
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